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Spectral imaging of the Central Molecular Zone in multiple 7-mm 
molecular lines 
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We have imaged 24 spectral lines in the Central Molecular Zone (CMZ) around the Galac- 
tic Centre, in the range 42 to 50 GHz. The lines include emission from the CS, CH3OH, 
HC 3 N, SiO, HNCO, HOCO+, NH2CHO, OCS+, CCS, C 34 S, 13 CS, 29 SiO, H 13 CCCN, 
HCC 13 CN and HC5N molecules, and three hydrogen recombination lines. The area covered 
is Galactic longitude —0.7 to 1.8 deg. and latitude —0.3 to 0.2 deg., including the bright cores 
around Sgr A, Sgr B2, Sgr C and Gl .6-0.025. This work used the 22-m Mopra radio telescope 
in Australia, obtaining ~ 1.8 km s _1 spectral and ~ 65 arcsec spatial resolution. We present 
peak images from this study and conduct a principal component analysis on the integrated 
emission from the brightest 10 lines, to study similarities and differences in the line distribu- 
tion. We examine the integrated line intensities and line ratios in selected apertures around the 
bright cores, as well as for the complete mapped region of the CMZ. We compare these 7-mm 
lines to the corresponding lines in the 3-mm band, for five molecules, to study the excitation. 
There is a variation in 3-mm to 7-mm line ratio across the CMZ, with relatively higher ratio 
in the centre around Sgr B2 and Sgr A. We find that the lines are sub-thermally excited, and 
from modelling with RADEX find that non-LTE conditions apply, with densities of order 10 4 
cm -3 . 

Key words: ISM:molecules - radio lines:ISM - ISM:kinematics and dynamics. 
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1 INTRODUCTION 

The Central Molecular Zone (CMZ) is the region within a few hun- 
dred parsecs of the Galactic Centre (lMorris & Serabvnl ll996) with 
dense molecular gas, of total molecular mass a few times 10 7 Mq 
jFerriere. Gillard & Jeanl 12007). The CMZ is an unusual environ- 
ment, with high molecular gas densities, large velocity dispersion 
and s hears, high temperature, and rich chemistry: see Uones et al.l 
d2012h and references therein. 

Although there is on-going star formation in the CMZ, the 
star-formation rate (SFR) is an order of magnitude less than 
would be expected fro m such a large amount of molecular gas 
(lLongmore et al.l l2012h using common SFR scaling relations. 
Hence, studying the properties of the molecular medium in the 
CMZ is important in understanding the enviromental dependence 
of star formation, particularly the environments of the centres of 
galaxies. 

We present here observations of the CMZ of 24 lines in the 
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7-mm band (including three hydrogen recombination lines), with 
resolution around 65 arcsec. 

This is part of a larger project in which have also imaged 
the CMZ area in multipl e lines in the 3-m m band, with reso- 
lution around 40 arcsec (IJones et al.l l2012h . and in the 12-mm 
band, with resolution around 120 arcsec, as pa rt of the wide r H2O 
southern Galactic Plane Survey (HOPS. Iwalsh et al.ll2008l 1201 ll : 
IPurcell etall 120121) . A further Mopra survey to image CO lines 
over a larger area of the CMZ is ongoing. We have also studied the 
smaller Sgr B2 sub-area of the CMZ with the molecular lines over 
the full 3-mm band from 82 to 1 14 GHz, and 7-mm band from 30 to 
50 GHz, in Jones et al. (2008a) and Uones et alJd201lh respectively. 

We use the distance to the Galactic Centre of 8.0 kpc 
dReid et aDl2009h . as per our previous papers, so 1 arcmin corre- 
sponds to 2.3 pc. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made with the 22-m Mopra radio tele- 
scope, of the Australia Telescope National Facility (ATNF) using 
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Table 1. The lines imaged here in the 42 to 50 GHz range. The rest frequen- 
cies with a * indicate the frequency used for lines corresponding to multiple 
transitions, some also indicated as group (gp) in the transition list. The mul- 
tiple components are, in general, blended due the large velocity widths in 
the CMZ area. The Mopra observations covered the whole 42 to 50 GHz 
range, so include wea ker lines than are listed here, particularly in Sgr B2, 



range, so include weake 
see l Jones etakl J201lh . 



Rough 

frequency 

GHz 



Line ID 
molecule 
or atom 



Transition 



Exact 

rest frequency 

GHz 



42.39 


NH 2 CHO 


2(0,2)- 1(0,1) gp 


42.386070* 


42.60 


HC 5 N 


16-15 


42.602153 


42.67 


HCS+ 


1-0 


42.674197 


42.77 


HOCO+ 


2(0,2) -1(0,1) 


42.7661975 


42.82 


SiO 


1 - v=2 


42.820582 


42.88 


29 SiO 


1 - v=0 


42.879922 


42.95 


H 


53 a 


42.951968 


43.12 


SiO 


1 - v=l 


43.122079 


43.42 


SiO 


1 - v=0 


43.423864 


43.96 


HNCO 


2(0,2)- 1(0,1) gp 


43.962998* 


44.07 


CH3OH 


7(0,7)- 6(1,6) A++ 


44.069476 


44.08 


H 13 CCCN 


5-4 


44.084172 


45.26 


HC 5 N 


17-16 


45.264720 


45.30 


HC 13 CCN 


5-4 


45.297346 




HCC 13 CN 


5-4 


45.301707* 


45.38 


CCS 


4,3 - 3,2 


45.379029 


45.45 


H 


52 a 


45.453719 


45.49 


HC 3 N 


5-4gp 


45.490316* 


46.25 


13 cs 


1-0 


46.247580 


47.93 


HC5N 


18-17 


47.927275 


48.15 


H 


51a 


48.153597 


48.21 


c 34 s 


1-0 


48.206946 


48.37 


CH30H 


1(0,1)- 0(0,0) A++ 


48.372467* 




CH30H 


1(0,1)- 0(0,0) E 


48.376889 


48.65 


ocs 


4-3 


48.6516043 


48.99 


cs 


1-0 


48.990957 



the UNSW MOPS digital filterbank. The Mopra MMIC (Mono- 
lithic Microwave Integrated Circuit) receiver has a bandwidth of 
8 GHz, and the MOPS backend can cover the full 8-GHz range 
simultaneously in the broad band mode. This gives four 2.2-GHz 
sub-bands each with 8192 channels of 0.27 MHz. The lines in the 
CMZ are broad, so that these channels, corresponding to around 
1.8 km s _1 , are quite adequate. 

The Mopra receiver covers the range 30 to 50 GHz in the 7- 
mm band. We chose the tuning centred at 45.8 GHz, to cover the 
range 41.8 to 49.8 GHz, giving the spectral lines summarised in 
Table Q] This range was chosen to include the strong lines of CS 
near 49.0 GHz and SiO near 43.4 GHz, and several other bright 
lines as listed in the Table. This is also the densest part of the 30 
to 50 GHz band for bright molecular lines, as seen by Ijones et all 
bOllh in Sgr B2. 

The area was observed wit h on-the-fly (OTF) mapping, in a 
similar way to that described in I Jones et al.1 (12012), for the 3-mm 
CMZ observations. 

The 10 x 10 arcmin 2 OTF maps used a scan rate of 8 arc- 
sec/second, with 24 arcsec line spacing, taking a little over an 
hour per block. We made pointing observations of SiO maser po- 
sitions (AH Sco or VX Sgr), after every second map, to correct 
the pointing to within about 10 arcsec accuracy. The system tem- 
perature was calibrated with a continuous noise diode, but the 7- 
mm calibration system (unlike the 3-mm system) does not include 



ambient-temperature paddle scans, so some extra steps in process- 
ing were needed (see below) to correct for the atmospheric opacity. 
We included observations of line sources (Sgr B2, M17SW and 
G345.5+1.5) a few times each day, under different weather condi- 
tions and at different elevations, to check the calibration. 

We used position switching for bandpass calibration with 
the off- source reference position observed before each 10 ar- 
cmin long source scan. The reference position ((a, 5) 32000 = 
17 h 51 m 03. s 6, -28°22 / 47 // , or I = 1.093 deg., b = -0.735 deg.) 
was carefully selected to be relatively free of emission. 

Each 10x10 arcmin 2 block was observed twice, with Galactic 
latitude and longitude scan directions, to reduce scanning direction 
stripes and to improve the signal to noise. The overall 2.5 x 0.5 
deg 2 area was covered with a 15 x 3 grid of blocks, separated by 
10 arcmin steps, making 90 OTF observations required. In practice 
a few more OTF observations were used, including some areas with 
OTF maps stopped by weather or other problems and re-observed, 
or just areas re-observed in better weather. The area covered is be- 
tween -0.72 to 1.80 degrees Galactic longitude, and -0.30 to 0.22 
degrees Galactic latitude. 

The observations were spread over two periods, in March 
2010 and 2011. The conditions, in the southern hemisphere au- 
tumn, are warmer and wetter than the best winter observing pe- 
riods. Although the blocks worst affected by poor weather were 
re-observed, the range of conditions during the different observa- 
tions means that there are inevitably some blocks with greater noise 
level than others. In practice, we observed the central area includ- 
ing Sgr B2 and Sgr A in 2010, and the outer areas in 2011. The 
weather in 2011 Mar was not as good as in 2010, so that the sys- 
tem temperature was higher, and the noise level greater in the outer 
areas. 

The OTF data were turned into FITS data cubes wit h the LIVE- 
DATA and GRIDZILLA packageo in a similar way to Ijones et all 
J2012I) . The bandpass correction used the off- source spectra, and a 
second order polynomial, and the gridding used a median filter. 

The FITS cubes were then read into the MIRIAD package for 
further processing and analysis. The velocity pixels were 0.27 MHz 
or ~ 1.8 km s _1 , but to improve the signal-to-noise, and ensure the 
lines were Nyquist sampled, we applied a 3 -point Hanning smooth- 
ing, leading to an effective resolution 0.54 MHz or ^3.6 km s _1 . 
This is quite adequate spectral resoution for the broad lines in the 
CMZ, which are > 10 km s _1 wide, but for narrow spectral fea- 
tures, such as the maser lines, we use the data cubes in the original 
unsmoothed versions. 

Small residual baseline offsets in the spectra were fitted and 
removed with a first order polynomial using the MIRIAD task con- 
tsub. 

The resolution of the Mopra beam varies between 1.14 ± 
0.05 arcmin at 43 GHz to 0.99 ± 0.04 arcmin at 49 GHz 
dUrquhart et al .1120101) . We take this as the beamsize inversely pro- 
portional to frequency, with 1.06 arcmin or 64 arcsec at 46 GHz in 
the centre of the frequency range observed. 

The main beam efficiency of Mopra varies between 0.53 =b 
0.01 at 43 GHz and 0.43 ± 0.01 at 49 GHz, and the extended 
beam effici ency between 69 ± 0.01 at 43 GHz and 0.56 ± 0.01 
at 49 GHz dUrq uhart et al .1 120 1 Oh . Since we are largely concerned 
in this paper with the spatial and velocity structure, we have mostly 
left the intensities in the T% scale, without correction for the beam 
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efficiency onto the Tmb or Txb scales, except for some of the 
quantitative analysis sections. 

A flux correction has been applied to the data cubes to 
take into account the atmospheric attenuation. The test spectra 
(of Sgr B2, M17SW and G345.5+1.5) showed a significant (a 
few times 10 %) decrease in flux associated with higher sys- 
tem temperatures. This was consistent with the effect expected 
with increased optical depth r in the atmosphere, at lower eleva- 
tion or higher humidity, increasing T sys with the addition of term 
Tamb(l - exp(-r)) rsj rT amb (where T amb ~ 290 K is the atmo- 
spheric temperature) and decreasing the flux by factor exp(— r). 
The system temperatures were measured in each of the 2 GHz sub- 
bands, across the overall 8 GHz range observed. When combined 
in the data cubes, averaged over the two passes in latitude and lon- 
gitude scanning directions, the means were 84, 93, 106 and 126 K 
at around 43, 45, 47 and 49 GHz respectively, with RMS scatter 
10, 11, 13 and 15 K respectively. However, due to outliers of data 
taken in poor weather, the full ranges were 68 - 190, 75 - 201, 85 
- 219 and 101 - 231 K respectively. This provides a measure of the 
average optical depth r as a function of position in the images, so 
allowing a correction for it to be applied. The flux correction fac- 
tors applied to the cubes, as a function of position derived from the 
T sys images were 1.150 ± 0.044, 1.180 ± 0.054, 1.233 ± 0.068 
and 1.276 ± 0.087 in the four sub-bands. 



3 RESULTS 

3.1 Line data cubes 

The frequencies and line identifications for the 24 strongest lines 
in the 41.8 to 49.8 GHz range, are listed in Table Q] These data 
are taken from the NIST online cataloguq 2 ] o f lines known in the 
interstellar medium (ILovas & Dragosetl 120041) and the splatalogue 
compilation^. 

Some statistics of the data cubes are listed in Table |2] The 
RMS noise ranges from 24 to 76 mK, increasing with frequency 
across the 8 GHz range, with the highest noise for the 47.93 GHz 
HC5N line, near the 47.8 GHz edge of the sub-band. 

The velocity range listed in Table indicates the range of sig- 
nificant emission (roughly greater than 3 a) detected in the data 
cubes. The line widths in the CMZ are large, and there is a large 
velocity gradient across the CMZ (in the deep potential well of the 
Galactic Centre). Hence the integrated emission needs to summed 
over a large velocity range. This makes the integrated emission 
quite sensitive to low level baselevel offsets which are small in 
terms of brightness (e.g. of order 10 mK) but become significant 
(e.g. of order K km s _1 ) when integrated over the velocity range 
(order 100 km s _1 ). The 7-mm Mopra data here have better base- 
level stability than the 3 -mm data of Uones et al.1 (120121) . due to the 
lesser sensitivity to the weather at lower frequencies. However, the 
baseline problems do provide a limit, along with the thermal noise, 
to the accuracy of the data. 

These data cubes are publicly available in the CSIRO-ATNF 
archive, accessible from http://atoa.atnf.csiro.au/CMZ together 
with the 3 -mm and 12-mm data from the CMZ and Sgr B2|j 



2 http://physics.nist.gov/PhysRefData/Micro/Html/contents.html 

3 http://www.splatalogue.net/ 

4 Further details may be also be obtained from the UNSW web site at 
w w w.phy s . uns w. edu . au/mopracmz . 



3.2 Peak emission images 

We show the 2-dimensional distribution of the line emission in Figs 
Q]to|4] We use peak images as these are are much more robust to the 
effect of low level baseline problems than the integrated emission 
images. They also better show fine spatial structure as they are less 
affected by multiple velocity components along the line of sight. 
The colour-scale or grey-scale displays start at the 3 a brightness 
level, so do not show noisy features in areas without significant 
line emission. However, low-level emission outside the areas high- 
lighted in these plots, may be found to be significant when the data 
are integrated over larger areas than the beam, or smoothed in ve- 
locity. 

There are some artifacts in some of these peak brightness im- 
ages, notably areas at longitude 359.6 deg. with higher noise, giving 
maxima above the overall 3 a level. 

The lines of CS, CH3OH (at 48.37 GHz), HC 3 N, SiO and 
HNCO are relatively strong and widespread (Fig. [T). The lines of 
HOCO+ (Fig. ID, NH2CHO, OCS, HCS+, and CCS (Fig.0 are 
weaker, and detected mostly in the Sgr B2 and Sgr A areas. The 
lines of HNCO, HOCO+ and OCS highlight the core of Sgr B2, 
and also the area of Gl. 6-0.025 (as discussed below in section [33l >. 

The weaker isotopologues of CS, SiO and HC3N, namely 
C 34 S and 13 CS in Fig. 29 SiO, H 13 CCCN and HCC 13 CN 
(blended with HC 13 CCN) in Fig. [3] are also detected around the 
Sgr B2 and Sgr A areas. Note that the line of H 13 CCCN at 44.08 
GHz (Fig. [3]) is highly confused with the line of CH3OH at 44.07 
GHz (Fig. |4]), so the peak image has additional strong CH3OH 
maser peaks which overlap the frequency range. In principle, these 
isotopologues could be used to measure isotope ratios (C/ 13 C, S/ 34 
and Si/ 29 ) near the Galactic Centre. However, in practice, the low 
signal noise of the mapping data for these weaker isotopolgues 
means that they are only detected at the strong cores, where the 
strong isotopologues (CS, SiO and HC3N) are optically thick, so 
the line ratios do not give the intrinsic isotope ratios. 

The three lines of HC5N (Fig. [3]) show emission from Sgr B2, 
but the higher noise at the higher frequency means that the weaker 
emission at Sgr A and Gl. 6-0.025 detectable in the two lower fre- 
quency lines is lost in the noise for the last line. 

The line of CH3OH at 44.07 GHz (Fig© has maser as well as 
thermal emission, so shows diffuse thermal emission around Sgr B2 
and Sgr A, plus a number of maser spots distributed over a wider 
area of the CMZ. 

The two maser lines of SiO in Fig.|4]are dominated by a strong 
maser peak (at I = 0.550 deg., b = -0.057 deg.) but also show 
a number of weaker maser spots. These are vibrationally excited 
lines, unlike the thermal v=0 line of SiO at 43.42 GHz, and so have 
a very different distribution to the latter. 

The three hydrogen a radio recombination lines are at Sgr B2, 
which is where the thermal bremsstrahlung emission in the CMZ is 
strongest. 

More details on the 7-mm line distributions in the area around 
Sgr B2 are given in ljones et al.l d201lh . 



3.3 Principal Component Analysis 

We can identify and quantify similarities and differences between 
the CMZ line emission here with Principal Component Analysis 
(PC A). This technique describes multi-dimensional data sets by lin- 
ear combinations of the data that describe the largest variance (the 
most significant common feature) and successively smaller vari- 
ances (the next most significant features). For integrated emission 
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Figure 1. Peak brightness images for the lines of CS, CH3OH, HC3N, SiO, HNCO and HOCO+. These are the six lines we observed with the most widespread 
strong distribution. The grey-scale is peak brightness as T^ in K here and in Figs[21to|4l 
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Figure 2. Peak brightness images for the lines of NH 2 CHO, OCS, HCS+ , CCS and the CS isotopologues C 34 S and 13 CS. 
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Figure 3. Peak brightness images for the lines of SiO isotopologue 29 SiO, HC 3 N isotopologues H 13 CCCN and HCC 13 CN blended with HC 13 CCN, and 
three transitions of HC 5 N. Note that the 44.08 GHz H 13 CCCN line is confused with the 44.07 GHz CH 3 OH maser line (Fig.gJ. 
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Figure 4. Peak brightness images for lines of CH3OH in a maser transition, two maser transitions of SiO, and three hydrogen a recombination lines. 
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Table 2. Statistics of the data cubes. Both the RMS noise and peak brightness temperature are in TJ, after the calibration correction for atmospheric attenuation. 
The velocity range is that with line emission significantly above the noise level. We also quote the position and velocity of the peak pixel. 



Line 


Molecule 


RMS 


Velocity 


Peak 


Peak 


Peak 


Freq. 


or atom ID 


noise 


range 




lat. long. 


vel. 


GHz 




mK 


km s — 1 


K 


deg. deg. 


km s _1 



42.39 


NH 2 CHO 


24 


-5,110 


0.30 


359.861 


-0.082 


8 


42.60 


HC 5 N 


25 


0,100 


0.33 


0.687 


-0.020 


65 


42.67 


HCS+ 


26 


-5,85 


0.24 


0.675 


-0.020 


71 


42.77 


HOCO+ 


27 


-5, 100 


0.39 


0.694 


-0.020 


62 


42.82 


SiO 


26 


-150, 200 


3.22 


0.550 


-0.057 


-25 


42.88 


29 SiO 


27 


0,85 


0.17 


0.836 


-0.248 


49 


42.95 


H53a 


29 


35,95 


0.39 


0.669 


-0.037 


64 


43.12 


SiO 


28 


-150, 200 


2.74 


0.550 


-0.057 


-26 


43.42 


SiO 


40 


-140, 205 


1.48 


359.893 


-0.076 


15 


43.96 


HNCO 


38 


-160, 185 


2.68 


1.649 


-0.053 


51 


44.07 


CH3OH 


30 


-60, 195 


3.37 


359.986 


-0.069 


47 


44.08 


H 13 CCCN 


30 


0,90 


0.22 


0.698 


-0.028 


61 


45.26 


HC 5 N 


30 


0,80 


0.32 


0.690 


-0.026 


63 


45.30 


HCC 13 CN 


29 


30, 105 


0.25 


0.675 


-0.020 


68 


45.38 


CCS 


31 


-5, 100 


0.35 


0.687 


-0.020 


70 


45.45 


H52a 


38 


30,90 


0.42 


0.669 


-0.037 


63 


45.49 


HC 3 N 


37 


-150, 200 


2.67 


359.981 


-0.078 


16 


46.25 


13 CS 


34 


-15, 110 


0.49 


359.983 


-0.071 


52 


47.93 


HC 5 N 


73 


20,80 


0.37 


0.675 


-0.022 


67 


48.15 


H51a 


67 


40,85 


0.45 


0.672 


-0.035 


58 


48.21 


C 34 S 


59 


-25, 140 


0.81 


359.982 


-0.073 


51 


48.37 


CH3OH 


76 


-160, 205 


3.76 


0.689 


-0.024 


71 


48.65 


ocs 


55 


-10, 105 


0.72 


0.691 


-0.029 


64 


48.99 


cs 


75 


-185,210 


3.51 


359.984 


-0.069 


54 



images here, PCA gives a series of images which contain the most 
significant features (Fig. [5} and a set of coefficients (Fig. [6) which 
describe how the individual line images are made up of linear com- 
binations of these PCA images. 

We have used this PCA technique on the Mopra line surveys of 
the C MZ at 3-mm (Jones et"aHl2012h . the Sgr B2 area at 3-mm and 



7-mm (Jones et al. 2008b 



and the G333 molecular cloud 



complex at 3-mm dLo et al.l2009h . These papers, and the references 



therein, give more details on the technique. 

Since the PCA uses normalised versions of the input data sets, 
it does not work well with low signal to noise data, as the noisy 
data is scaled up, causing the result to be dominated by spurious 
features. Hence we use the ten strongest lines here (CS, CH3OH 
(at 48.37 GHz), HC 3 N, SiO, HNCO, HOCO+, NH 2 CHO, OCS, 
C 34 S and 13 CS) out of the total of 24 lines. We consider the whole 
observed area, less a small cropping at the low sensitivity edges of 
the are a. This is different to the PCA analysis of the CMZ 3-mm 
lines in ljones et al.1 d2012h , where we applied an additional mask of 
low surface brightness areas, leaving only the stronger areas around 
Sgr B2 and Sgr A, because the 3-mm data were more badly affected 
by baseline ripples. 

We use the versions of the integrated line emission (over the 
velocity ranges given in Table |2} with clipping below the 3 a level, 
as these are less affected by baseline ripple than the integrated emis- 
sion without clipping. However, we have checked the PCA results 
from integrated emission, with and without this clipping, and find 
qualitatively similar results. Hence the PCA results are not too se- 
riously affected by the small amount of missing flux caused by the 
3 cr-level clipping. 

The images of the three most significant components of the 
integrated emission are shown in Fig. [5] These three components 



describe 68.3, 11.5 and 7.3 percent respectively of the variance in 
the data. The normalised projection of the molecules onto the prin- 
cipal components are shown in Fig. [6] with successive pairs of com- 
ponents. 

The first principal component (Fig. [5} shows emission com- 
mon to all ten lines. As all of these lines are quite similar in overall 
morphology, it describes a large fraction (68.3 percent) of the vari- 
ance. This behaviour is similar to that seen at 3-mm (jjones et al.l 
120121) . where the first component displays an averaged intensity 
image of all the molecular lines (albeit with a factor two better res- 
olution). 

The second principal component (Figs. [5] and [6) shows the 
major difference (1 1.5 percent of the variance) among the ten lines, 
which is the difference between the brightest regions (blue in the 
colour version of Fig. [5] dark in the black and white version) and the 
fainter emission (red and yellow in the colour version, light in the 
black and white version). We attribute this to the dense cores hav- 
ing high optical depth in CS, SiO and (to a lesser extent) CH3OH, 
HNCO and HC3N, so that the peaks have less emission than would 
be expected from optically thin lines (as the others are). This is 
borne out by comparing specifically CS to the isotopologues C 34 S 
and 13 CS. This pattern of the second principal compo nent show- 
ing the optical depth effect is similar to that found in Jones et al. 
(2012) for the 3-mm CMZ lines. It is the case that this second prin- 
cipal component is related to the strength of the line, so there could 
be some concern about it being affected by the missing flux in using 
the integrated emission with 3 a clipping, however (as noted above) 
the PCA using the integrated emission without clipping gives sim- 
ilar results. 

The third principal component shows the next most signifi- 
cant difference among the ten lines (7.3 percent of the variance). 



© 2013 RAS, MNRAS 000,[Q{T5] 



7-mm spectral imaging of the CMZ 9 



This component shows differences between the lines of CS, C 34 S 
and 13 CS on one side, with a relative excess around Sgr A (red 
and white in the colour version of Fig. [5] dark in the black and 
white version) and the lines of HOCO+, HNCO and OCS on the 
other, which show a relative excess around Gl. 6-0.025 and the 
peak of Sgr B 2 (blue in the colour version, light in the black and 
white version). The relative excess of HOCO + , HNCO and OCS 
around Gl. 6-0.025 can be seen by eye in Figs. Q] and |2] (albeit 
plotted as peak emission, not integrated emission). HNCO has in- 
deed been shown to vary significantly in r elative intensity to d ense 
gas tracers such as CS across the CMZ b y Martin et al. (2008|) and 
lAmo-Baladron. Martin-Pintado. & Martini (1201 ll) . with the sugges- 
tion that this provides a diagnostic tool to distinguish between re- 
gions dominated by shock or by PDR emission from the molecules. 
In the 3 -mm PC A analysis (IJones et al.ll20121) the shock- sensitive 
species HNCO and SiO were distinguished from the dense gas trac- 
ers. 

The line differenc es in the Sgr B2 are a at 7-mm and 3 -mm ar e 
shown in more detail in ljones et al] (1201 11) and Uones et al.l(12008bl) . 
and discussed there, but HOCO + and HNCO do highlight a chem- 
ically distinct region in the Sgr B2 area. 

We note that after the "optical depth correction" effect of 
the second principal component that CS and weaker isotopologues 
C 34 S and 13 CS now show similar behaviour in the third principal 
component (same sign, 13 CS and CS similar magnitude), giving 
more confidence that this third component is showing real differ- 
ences in the molecular distributions. 

We also note that the PCA images from the 7-mm lines here 
(Fig. [5]) for the first three components agree very well w i th the PCA 
images from the 3 -mm lines in figure 9 of IJones et al.l (J2012I) . de- 
spite the fact that the latter are derived from different lines. This 
gives confidence that the PCA technique is identifying physically 
meaningful features. The sign of the principal components i s arbi- 
trary, and the grey- scale display of figure 9 of Uones et al.l (120121) 
was chosen differently to Fig.[5]here, so the visual displays of the 
two sets of PCA images are a bit different, but the data are quite 
similar. 



3.4 Line intensities and line ratios 



Table 3. Apertures selected for analysis. The areas are rectangles in Galactic 
coordinates, centred on I and b. The projected area on the sky assumes a 
distance of 8.0 kpc. 



The previous subsection (13.31 ) indicates that there are differences in 
the ratio of lines for different areas within the CMZ. We can quan- 
tify this by integrating the line emission over different apertures. 
The apertures chosen are listed in Tableland plotted in Fig. [7] The 
six regions are the entire region of the CMZ that we have mapped, 
the four main emission features namely Sgr A, Sgr B2, Sgr C and 
G1.3, and a smaller region around the core of Sgr B2 . They are 
effectively the same as those used in I Jones et al.1 d2012h for similar 
analysis of the CMZ in 3-mm lines, with very minor changes due to 
the data at 7-mm being on different (coarser) pixels to the previous 
3-mm data. 

Table HI presents the integrated line fluxes for these apertures, 
as ff Tb dv &A in K km s _1 pc 2 . The brightness temperature Tb is 
the extended beam te mperature, corrected for the extended aperture 
efficiency of ~ 0.65 ( Urquhart et al .1120101) . as discussed below. 

The line emission is integrated over velocity, with the accuracy 
limited by low level baseline ripples in the spectra, as discussed in 
subsection [TT] To reduce this problem, we integrated over a veloc- 
ity range which covers the expected emission while excluding too 
much extra data in the cubes. Hence this velocity range varies with 
the area (due to differential velocity across the CMZ, and different 
velocity widths in the different areas) and the line. Uncertainties 



Region 


I 


b 


Width 


Height 


Area 


Area 




deg 


deg 


arcmin 


arcmin 


arcmin 2 


pc z 


CMZ 


0.545 


-0.035 


151.2 


30.0 


4536 


24 600 


Sgr A 


0.042 


-0.055 


30.8 


13.2 


406 


2200 


SgrB2 


0.735 


-0.045 


29.2 


18.4 


537 


2910 


SgrC 


-0.462 


-0.145 


20.0 


16.8 


336 


1820 


G1.3 


1.325 


0.015 


24.0 


20.8 


499 


2700 


SgrB2 


0.675 


-0.025 


6.0 


6.4 


38.4 


208 


Core 















are quoted, based on observed integrated intensities fluctuations, 
caused by these baseline offsets, over parts of the data cubes with- 
out lines. However, we do note that the baseline ripples vary be- 
tween the different lines, as does the noise in the cubes (Table |2}, 
so that these uncertainty estimates are only indicative. 

These line luminositie s can be compared to those of 3-mm 
lines over the same areas in Uones et al.l (120121) . and are quoted in 
the form of ff Tb dv dA for ease of comparison with other galax- 
ies. 

We do not include all the lines in this table, but exclude the 
weakest lines which are generally below the 3 a level of sig- 
nificance in the areas given the baseline problems (H 13 CCCN, 
HCC 13 CN, HC 5 N, SiO masers, RRLs). 

Table \5\ presents the integrated line fluxes over the six aper- 
tures, normalised by the value for the CS line. The Sgr B2 area 
generally has an enhanced line flux ratio to CS, as compared to the 
other areas, and the sub-area Sgr B2 Core has even more enhanced 
line ratios. While we expect Sgr B2 to be chemically rich, the high 
ratios relative to CS are likely to be in part due to the high optical 
depth of CS in this area. This can be checked with the isotopologues 
13 CS and C 34 S. These do indeed show higher C 34 S/CS ratio in Ta- 
ble [5] at Sgr B2, and particularly at the Core, which indicates high 
optical depth in CS. Expressed in the more commonly used form 
of the ratio CS/C 34 S (i.e. the inverse to that tabulated) this is lower 
ratio CS/C 34 S in Sgr B2 (8.8) and the Sgr B2 Core (7.5) compared 
to the other regions. The situation in the weaker 13 CS line is less 
clear, with lower CS/ 13 CS at the Sgr B2 Core, of 14 compared to 
the expected optically thi n value of 24 for the Galactic Centre from 
lLanger & Penziasl dl990L but overall ratio 20 in the larger Sgr B2 
region (and an anomalously high ratio 63 for Sgr C, which may 
indicate baselevel problems). 



3.5 Ratio of 7-mm and 3-mm lines 

The observed integrated line intensity f Tdv of a spectral line can 
be related to the column density N u of molecules in the upper level 
of the transition, for optic ally thin lines using standard as sumptions 
of LTE radiative transfer (IWilson. Rohlfs & Huttemeisterl2009l) . by 
N u = (87rv 2 k J he 3 Aui) f Tsdv, where A u \ is the Einstein coef- 
ficient. The total column density N of the molecule is given by 
N = (N u /g u )Q T exp(E u /kT ex ) where QT(T ex ) is the partition 
function, T ex is the excitation temperature, E u is the energy of the 
upper level and g u is the statistical weight of the upper level. 

The excitation temperature T ex describes the relative popula- 
tions in different levels, and is typically not the same as the gas 
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Figure 5. The first three principal component images derived from the ten strongest lines (CS, CH 3 OH (at 48.37 GHz), HC3N, SiO, HNCO, HOCO+, 
NH2CHO, OCS, C 34 S and 13 CS), which describe 68.3, 11.5 and 7.3 percent of the variance respectively. The first principal component describes the common 
features of the ten lines, and the second principal component describes the most significant differences between the lines. The further principal components 
describe successively smaller differences. 
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Figure 6. The component vectors quantifying how the integrated line images are composed of the sum of different scalar amounts of the of the principal 
component images (Fig.|5j- 



kinetic temperature Tkin, as at low d ensities the excitatio n depends 
on the radiation as well as collisions JWilson et all 2009). 



As we have 3 -mm line observations of the CMZ from 
I Jones et al] d2012h to match the 7-mm observations here, we can 
consider the ratio of integrated line intensities f Tdv, which is re- 



lated (subject to assumptions, and caveats below) to the excitation 
temperature T ex - 

We consider five molecules, that is HNCO, HOCO+, HC 3 N, 
SiO and 13 CS, common to both our 7-mm and 3 -mm CMZ spectral 
imaging. The 3-mm data were smoothed to the resolution of the 7- 
mm data, using the beamsizes (62 to 68 arcsec) interpolated from 
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Figure 7. The six areas selected for analysis: around Sgr A, Sgr B2, Sgr C, G1.3, a smaller area at the Sgr B2 core, and the larger CMZ area, plotted on the 
peak CS emission. 

Table 4. Integrated line fluxes as ff Tb dv <±A, integrated over velocity, and the areas given in Table [3] The integrated fluxes are in units 10 3 K km s _1 pc 2 , 
The 1 a uncertainties are given, and the upper limits are at the 3 a level. 



Line 



CMZ 



Sgr A 



SgrB2 



SgrC 



G1.3 



Sgr B2 Core 



CS 


1590 ± 80 


218 ±11 


297 ± 15 


157 ±8 


246 ± 12 


33.2 ±1.7 


CH 3 OH 


650 ± 30 


66 ±4 


182 ±9 


48 ± 2.5 


130 ±7 


29.4 ± 1.5 


HC3N 


287 ± 15 


33.2 ± 2.2 


82 ±5 


16.6 ± 1.0 


57 ±3 


15.1 ±0.8 


SiO 


283 ± 15 


24.8 ± 1.9 


71 ±4 


14.4 ± 1.0 


72 ±4 


8.5 ± 0.5 


HNCO 


312 ± 16 


23.2 ± 1.8 


80 ±5 


14.6 ± 1.0 


66 ±3 


12.9 ± 0.7 


HOCO+ 


30 ±5 


<4 


8.0 ± 2.5 


<1.9 


5.5 ± 0.9 


1.69 ±0.28 


NH2CHO 


20 ±5 


<4 


7.9 ± 2.5 


<1.9 


4.1 ±0.9 


1.36 ±0.28 


OCS 


60 ±6 


<4 


25.2 ± 2.7 


<1.9 


16.6 ± 1.2 


3.6 ± 0.3 


HCS+ 


27 ±5 


5.2 ±1.4 


<7 


<1.9 


3.9 ± 0.9 


1.11 ±0.28 


CCS 


34 ±6 


<4 


13.5 ±2.5 


<1.9 


8.4 ± 0.9 


2.00 ± 0.29 


c 34 s 


114 ±8 


9.7 ±1.5 


33.5 ± 2.9 


10.5 ± 0.8 


19.8 ± 1.3 


4.4 ± 0.4 


13 cs 


56 ±6 


11.6 ± 1.5 


14.8 ± 2.5 


2.5 ± 0.6 


9.3 ±1.0 


2.4 ± 0.3 


29 SiO 


22 ±5 


<4 


<7 


2.0 ± 0.6 


3.0 ± 0.9 


<0.8 


CH3OH maser 


53 ±6 


9.6 ± 1,5 


12.0 ± 2.5 


4.0 ± 0.6 


10.1 ± 1.0 


4.0 ± 0.3 

















the Mopra 7-mm measurements of lUrquhart et al.l (120101) . The 3 
mm and 7-mm data expressed as T\ were converted to brightness 
temperatures Tb appropriate for extended emission which fills both 
the main beams and inner sidelobes, by dividing by the extended 
beam efficiencies at 3-mm (0 . 62 to 0.65 ) and 7-mm (0.66 to 0.69) 
interpolated from lLadd et all h005) and Urauhart et al.l (1201 Oh re- 
spectively. 

The ratio image of f Tsdv for the 3-mm and 7-mm lines was 
taken by summing 3 -dimensional pixels (voxels) over velocity in 
both cubes, over the region where both lines were above the respec- 
tive 3-cr level, and then dividing the integrated emission images. 
This reduces the effect of including data without significant emis- 
sion which adds artifacts due to low-level baseline offsets (particu- 
larly in the 3-mm data which is more weather affected), but should 
not bias the ratio due to missing flux, as we consider the ratio of 
emission from the same voxels. It is effectively a weighted mean 
of the ratio, considered on a voxel basis, with the weighting set to 
zero below the 3-cr level. 

The ratio images are presented in Fig. [8] as the nat- 
ural logarithm of the ratio of 3-mm to 7-mm fTsdv, 
HI T B ,3mmdv/ f T B ,7mmdv). The HNCO, HOCO+ and HC 3 N 
images show variations in the ratio across the CMZ, in the sense 
that the 3-mm to 7-mm ratio is higher around Sgr B2 and Sgr A, 
and is less around longitudes 1.0 deg. to 1.8 deg. There are some 
signs of similar variation in the ratio in the SiO image, but less 
clearly. 

The ratios are plotted in Fig. [9] averaged over latitude, to show 



the longitude variations more clearly. The plots for HOCO + and 
13 CS are noisy for longitude ranges with little data, but the HOCO + 
does show the difference in ratio between Sgr B2 and higher lon- 
gitudes. The plots for HNCO, HC3N and SiO do show ratio varia- 
tions. 

At the centre of Sgr B2 in Fig. [8] there is a small region of 
higher 3-mm to 7-mm ratio in HC3N, SiO and 13 CS, associated 
with the free-free radio peak of Sgr B2. We attribute this to the 
radiative transfer effect of line emission being absorbed in front 
of the strong continuum. This occurs to a greater exten t in the 7- 
mm band compared to the 3-mm band, as discussed in I Jones et al.l 
(120120 , leading to the enhanced 3-mm to 7-mm ratio. 

The line ratio differences can largely be interpreted as exci- 
tation temperature differences, as discussed above, with the caveat 
that other effects may be occuring, such as the continuum absorp- 
tion and signifcant optical depth (subsection l3.3i >. The line frequen- 
cies (v), Einstein coefficients (A u i), statistical weights (g u ) and up- 
per energy levels (E u ) were obtained from the splatalogue online 
compilation ( http://www.splatalogue.net/). 

The statistics of the line ratios 

\n(f TB,3mmdv/ f TB,7mmdv) (mean and standard devia- 
tion) are listed in Table [6] for the areas in Table |3^nd Fig. [7] The 
corresponding excitation temperatures T ex assuming optically thin 
LTE conditions are listed in Table |7] 

Note that these calculated excitation temperatures T ex are 
much lower than the kinetic temperature Tkin ~ 30 K. In practice, 
for the likely physical conditions in the CMZ, the energy levels 



© 2013 RAS, MNRAS 00Q,fTlfl5l 



12 P. A. Jones et al. 

Table 5. Integrated line fluxes, integrated over velocity and the areas given in Table [3] normalised to CS. The 1 a uncertainties are given, and the upper limits 
are at the 3 a level. 



Line 


CMZ 


SgrA 


SgrB2 


SgrC 


G1.3 


Sgr B2 Core 


CS 


1.00 def. 


1.00 def. 


1.00 def. 


1.00 def. 


1.00 def. 


1.00 def. 


CH 3 OH 


0.41 ± 0.03 


0.301 ± 0.022 


0.61 ± 0.05 


0.306 ± 0.022 


0.53 ± 0.04 


0.89 ± 0.06 


HC 3 N 


0.181 ±0.013 


0.152 ±0.013 


0.278 ± 0.021 


0.106 ±0.008 


0.230 ± 0.017 


0.46 ± 0.03 


SiO 


0.178 ±0.013 


0.113 ±0.010 


0.238 ± 0.019 


0.092 ± 0.008 


0.291 ± 0.021 


0.255 ± 0.020 


HNCO 


0.197 ±0.014 


0.106 ±0.010 


0.270 ± 0.021 


0.093 ± 0.008 


0.270 ± 0.019 


0.39 ± 0.03 


HOCO+ 


0.019 ± 0.004 


< 0.019 


0.027 ± 0.008 


< 0.012 


0.022 ± 0.004 


0.051 ±0.009 


NH2CHO 


0.013 ± 0.003 


< 0.019 


0.027 ± 0.008 


< 0.012 


0.017 ± 0.004 


0.041 ± 0.009 


OCS 


0.038 ± 0.004 


< 0.019 


0.085 ± 0.010 


< 0.012 


0.067 ± 0.006 


0.109 ±0.011 


HCS+ 


0.017 ± 0.004 


0.024 ± 0.007 


< 0.024 


< 0.012 


0.016 ± 0.004 


0.033 ± 0.008 


CCS 


0.021 ± 0.004 


< 0.019 


0.046 ± 0.009 


< 0.012 


0.034 ± 0.004 


0.060 ± 0.009 


c 34 s 


0.072 ± 0.006 


0.044 ± 0.007 


0.113 ±0.011 


0.067 ± 0.006 


0.080 ± 0.007 


0.133 ±0.012 


13 cs 


0.035 ± 0.004 


0.053 ± 0.007 


0.050 ± 0.009 


0.016 ± 0.004 


0.038 ± 0.004 


0.073 ± 0.010 


29 SiO 


0.014 ± 0.003 


< 0.019 


< 0.024 


0.013 ± 0.004 


0.012 ± 0.003 


< 0.024 


CH3OH maser 


0.033 ± 0.004 


0.044 ± 0.007 


0.041 ± 0.009 


0.025 ± 0.004 


0.041 ± 0.004 


0.121 ±0.012 



of these molecules are unlikely to be in LTE (see below, section [4]). 
Hence the values of T ex calculated here are a useful way to describe 
the relative populations of the levels that give rise to the 7-mm and 
3-mm lines, but should not be over-interpreted as applying to other 
levels. 

We note the line ratios are smaller beyond Sgr B2, and thus 
outside the 10 7 M rotating dust ring, that extends from Sgr B2 to 
Sgr A, citepmo+11. This is particularly noticeable with the G1.3 
dust core, where the higher excitation 3-mm lines are relatively 
weaker. 



4 DISCUSSION 

The kinetic temperat ure Tkin and hydr ogen density n in the CMZ 
have been fitted bv INagai et all (120071) using CO lines ( 12 CO 1- 
0, 12 CO 3-2, 13 CO 1-0) with a Large Velocity Gradient (LVG) 
model, giving typical values Tkin = 30 K and n = 3x 10 3 cm -3 . 

For such physical condition s Tkin and n we can con firm with 
the non-LTE model RADEX of Ivan der Tak et all d2007h . that the 
7-mm and 3-mm lines studied in section [331 give a low effective 
excitation temperature T ex - The hydrogen density n is much lower 
than the critical density n cr it at which the collisions would ther- 
malise the molecular levels to Tkin- We used the online version of 
RADEXQ to explore the parameter space that would reproduce the 
7-mm line fluxes and the 3-mm to 7-mm line ratios. We do note that 
this still assumes an isothermal homogeneous medium. 

The input parameters to RADEX are excitation conditions, 
Tkin, n and background temperature Tb g , and radiative transfer 
parameters, N and line width A v. If we fix Tkin — 30 K, 
A v = 10 km s _1 which is typical of these lines in the CMZ, and 
Tbg = 2.73 K for the microwave background temperature, then we 
have two free parameters n and N. Varying n changes the relative 
populations of the molecular levels, and hence the ratio of 3-mm to 
7-mm transitions. Varying N changes the line brightness tempera- 
ture Tr. With two free parameters, we can fit the two observable 
parameters, but only if the combination of observable parameters 
is in the range of the model, so this is not a trivial process. In prac- 
tice, varying one of n and N changes both the ratio and Tr due to 
non-linearities, such as optical depth. 

5 http://www.sron.rug.nl/vdtak/radex/radex.php 



The results of such modelling is given in Table [8] for four 
molecules. HOCO + is not one of the molecules with data files 
available for RADEX. We varied n and N to obtain typical val- 
ues of \n(T B ,3mm/T B ,7mm) and Tr (for the 7-mm line) given in 
the data in subsection [33] For each molecule, we present four fits, 
giving a typical Tr with mean ratio, higher and lower ratios by 0.3 
in the ln(ratio) or factor 1.35, and a fit for mean ratio and higher 
Tr. 

The derived values of n are a few times 10 4 cm -3 , which 
is somewha t larger than the t ypical value n — 3 x 10 3 cm -3 
obtained by INagai et all d2007h for CO, but these molecules with 
higher n cr it preferentially trace higher demsities. 

We have shown here that the observed ratio of the 3-mm to 
7-mm lines, and the line brightness temperatures, can be fitted with 
plausible values of n and N. With only two lines, we can only con- 
strain two parameters at a time, with other RADEX input parame- 
ters fixed. In practice, it is quite likely that the kinetic temperature 
Tkin and background temperature Tb g vary across the CMZ as well, 
so the actual variation in physical conditions is more complcated. 

The higher 3-mm to 7-mm line ratios around Sgr B2 and 
Sgr A, compared to longitudes around G1.3, can be partly ex- 
plained by higher density n, which means that the sub-thermal ex- 
citation gets closer to thermalisation with kinetic temperature Tkin- 
As noted above, however, from the data and analysis here we can- 
not rule out changes in Tb g and Tkin too. 



5 SUMMARY 

We have mapped a 2.5° x 0.5° region of the Central Molecular 
Zone (CMZ) of the Galaxy using the Mopra radio telescope in 
21 molecular and 3 hydrogen lines emitting from 42-50 GHz. The 
maps have spatial resolution ~ 65 arcsec and spectral resolution 
~ 1.8 km s -1 . For five species the emission is particularly strong 
and widespread; CS, CH3OH, HC 3 N, SiO and HNCO. For the 
other molecules emission is generally confined to the bright dust 
cores around Sgr A, Sgr B2 and G1.3. These data add to, and com- 
plement, previous studies we have undertaken in the 85-93 and 20- 
28 GHz portions of the mm- spectrum. As also seen in these other 
bands the molecular emission at 42-50 GHz is both widespread, 
but asymmetric, about the Galactic Centre, and turbulent. 

Principal component analysis (PCA) of the ten brightest emis- 
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Table 6. Mean and standard deviation of the ratio of 3-mm to 7-mm lines, expressed as ln(ratio), ie ln(T3 mm /T7 mm ) in the selected areas. 



Molecule 


CMZ 


SgrA 


SgrB2 


SgrC 


G1.3 


Sgr B2 Core 


HNCO 


0.07 ±0.32 


0.37 ±0.24 


0.43 ±0.21 


0.16 ±0.25 


-0.14 ±0.13 


0.58 ±0.11 


HOCO+ 


-0.04 ±0.33 


0.28 ±0.29 


0.07 ±0.29 


-0.12±0.15# 


-0.13 ± 0.14 # 


0.39 ±0.23 


HC 3 N 


-0.71 ±0.32 


-0.48 ±0.24 


-0.63 ±0.28 


-0.58 ±0.25 


-0.83 ±0.28 


-0.45 ±0.23 


SiO 


-0.30 ±0.18 


-0.22 ±0.14 


-0.25 ±0.15 


-0.21 ±0.21 


-0.34 ±0.14 


-0.25 ±0.13 


13 cs 


-0.36 ±0.25 


-0.30 ±0.16 


-0.37 ±0.17 


-0.23 ±0.18 


-0.35 ± 0.16 # 


-0.31 ±0.22 



# = small number of pixels (< 100) with good data (clipped at 3 a.) 



Table 7. The excitation temperatures T ex in K, assuming optically-thin LTE conditions, derived from the ratio of 3-mm to 7-mm lines. Note, however, that 
radex modelling in section [4] shows that assuming optically-thin LTE is too simplistic. 



Molecule 


CMZ 


SgrA 


SgrB2 


SgrC 


G1.3 


Sgr B2 Core 


HNCO 

HOCO+ 

HC 3 N 

SiO 

13 cs 


5.61H 

5-o±i-.g 

8-911:? 

9 c:+0.3 

9 9+0.4 
z - z -0.3 


7 o + 2.3 

6.6±?-| 

9.4111 

2-6±8:i 

9 9+O.2 
z - z -0.2 


7 7 +2.2 

'•'-1.4 

°-°-i.o 
8.71H 

2-6±8:| 

9 9+O.2 
z - z -0.2 


6.o±i:S 

4.8^+ # 

s.9±i-;g 

9 O+0.3 
z -°-0.2 


A Q+0.5 

'• y -0.9 
9 4+O.2 

9 9+O.4 # 
ZlZ -0.3 # 


7 9+2.2 
'• Z -1.4 

9.51H 
2.6±°;1 

9 9+0.3 
ZlZ -0.3 



# = small number of pixels (< 100) with good data (clipped at 3 a.) 



Table 8. Results of RADEX non-LTE models to reproduce typical 3-mm 
to 7-mm line ratios, and 7-mm line brightness observed in the CMZ. These 
are labelled m for mean or mid-range, 1 for low and h for high. We assume 
fixedT fcin = 30 K, Av = lOkms -1 and T bg = 2.7K. 



Molecule 


n 


N 


ln(T 3mm /T^mm ) 


Tr 




10 4 cm" 3 


10 14 cm" 2 




K 


HNCO 


0.65 


1.7 


-0.24 1 


1.0 m 


HNCO 


1.1 


1.7 


0.08 m 


1.0 m 


HNCO 


1.7 


1.8 


0.37 h 


1.0 m 


HNCO 


1.2 


4.5 


0.08 m 


2.5 h 


HC 3 N 


1.8 


0.50 


-1.02 1 


1.0 m 


HC 3 N 


2.7 


0.50 


-0.70 m 


1.0 m 


HC 3 N 


3.7 


0.50 


-0.41 h 


1.0 m 


HC 3 N 


2.5 


1.3 


-0.72 m 


2.4 h 


SiO 


1.5 


0.34 


-0.61 1 


0.49 m 


SiO 


4.2 


0.20 


-0.30 m 


0.52 m 


SiO 


8.0 


0.17 


-0.00 h 


0.51 m 


SiO 


4.5 


0.50 


-0.31 m 


1.2 h 


13 cs 


0.50 


0.30 


-0.66 1 


0.17 m 


13 cs 


1.5 


0.15 


-0.36 m 


0.17 m 


13 cs 


3.0 


0.13 


-0.06 h 


0.17 m 


13 cs 


1.5 


0.37 


-0.38 m 


0.40 h 



Line fluxes and luminosities are presented for the brightest 
dust cores, as well as for the entire CMZ, to aid in comparison with 
observations of nuclei of external galaxies where such components 
would be unresolved. 

Combining the data set with the higher-excitation lines we 
have mapped in the 3 mm band for five of the species (HNCO, 
HOCO + , HC 3 N, SiO and 13 CS) which indicates variations in the 
line ratios across the CMZ, with larger 3-mm to 7-mm line ratios 
observed outside the 10 7 M rotating dust ring, in particular for 
the molecular cloud complex G1.3. This implies cooler tempera- 
tures for the molecular gas here than in the dust ring. 

All these species are sub-thermally excited across the CMZ, 
with excitation temperatures considerably less than the kinetic tem- 
perature. Modelling of the emission using RADEX suggests gas 
densities for these molecules are, in general, a few x 10 4 cm -3 , a 
factor of a few higher than corresponding densities derived from 
CO line measurements. 

The full data set, comprising the data cubes for the 
21 molecular plus 3 hydrogen lines, is publicly available 
through the Australia Telescope National Facility archives 
(http://atoa.atnf.csiro.au/CMZ). A final data set in this series, con- 
sisting of data on the isotopologues of the CO J=l-0 lines over a 
5° x 1.5° of the CMZ, is currently being prepared for publication. 



sion lines quantifies the similarities and differences, with the first 
component, representing an averaged integrated map of all the 
molecular lines, contributing 68% of the variance shown in the 
data. The second component contributes 12% to the variance and 
highlights some differences in the dust cores, attributable to the 
brightest lines being optically thick in them. The third component 
contributes 7% to the variance and highlights some chemical dif- 
ferences, with HOCO + , HNCO and OCS being particularly promi- 
nent in G1.6 and Sgr B2. 
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Figure 8. The distribution of the ratio of integrated line brightness of 3-mm to 7-mm lines of HNCO, HOCO+, HC3N, SiO and 13 CS expressed as 
ln(J T B , 3mm dv/ f T B , 7mm dv). 
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Figure 9. The ratio of integrated line brightness of 3-mm to 7-mm lines of HNCO, HOCO+, HC3N, SiO and 13 CS expressed as 
In (J" TB,3mmdv/ f TB,7mmdv), averaged over latitude, to show the variation with longitude. Note the variation in line ratio across the CMZ, most clearly 
shown in HNCO, HC3N and SiO. 
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